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1.  INTRODUCTION 

In  recent  times  considerable  knowledge  has  been  gathered  that  relates  to 
both  the  methods  of  processing  signals  and  the  performance  that  results  when 
an  FM/cw  system  (radar)  operates  with  a  point  target.  The  simplest  systems 
use  a  single  down-conversion  receiver  (or  the  envelope  detector  equivalent 
receiver)  and  sense  the  resulting  Doppler  signal.  More  sophisticated  systems 
use  intermediate  frequency  (IF)  references  for  a  second  down-conversion  (or 
correlation)  operation;  these  references  can  be  either  simple  sinewaves  or 
"functional  references"  designed  to  maximize  system  response  for  an  anticipated 
specific  target  range.  These  correlation  operations  can  even  be  combined  with 
other  signal  processing  operations  to  produce  a  system  sensitive  to  the  direction 
of  the  Doppler  (incoming  versus  outgoing  target) .  Other  techniques  exist  for 
controlling  system  response  versus  target  range  and  even  reducing  undesirable 
range  sidelobes.  In  most  cases  the  analysis  of  these  advanced  systems  has 
been  developed  only  for  point  targets. 

Many  practical  applications  of  FM/cw  radars  call  for  the  system  to 
operate  above  the  earth's  surface.  In  fact,  the  surface  may  actually  be  the 
radar's  target.  For  example,  an  FM/cw  radar  functioning  as  a  fuzing  device 
may  be  required  to  sense  when  the  radar  has  reached  a  specific  height  above 
the  surface  so  as  to  detonate  its  charge.  In  this  case  the  surface  back- 
scatter  from  directly  below  the  radar  contains  the  most  important  information. 
However,  backscattered  power  from  other  directions  now  becomes  important 
because  the  surface  is  extended  and  is  not  a  point  target.  Analysis  of 
system  performance  becomes  more  difficult  with  extended  targets  and  the 
available  knowledge  base  is  not  as  well  developed  as  for  point  targets,  even 
for  simple  systems. 


In  this  report  we  examine  the  response  of  a  relatively  simple  form  of 
FM/cw  radar  when  moving  toward  a  homogeneous  flat  extended  (earth)  surface. 
The  assumption  of  homogeneity  does  not  prevent  us  from  considering  various 
types  of  surfaces,  such  as  grassy  fields,  water,  etc.;  it  only  prevents  us 
from  mixing  the  types . 

The  system  response  that  is  studied  is  the  output  (Doppler)  signal  from 
a  simple  down-conversion  system  that  uses  a  functional  reference.  The  power 
spectrum  of  the  response  is  found  and  the  form  of  the  reference  is  specified 
that  "optimizes"  the  power  spectrum. 

Results  found  herein  are  considered  preliminary  in  that  they  represent 
only  the  initial  findings  of  an  ongoing  study. 


In  the  development  of  our  principal  results  for  the  distributed  surface 
it  is  helpful  to  first  formulate  the  single,  or  point,  target  problem  in 
a  way  that  will  establish  a  particular  viewpoint. 
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V  V 


2.  POINT  TARGET  CASE 


Two  principal  purposes  are  served  by  Initially  formulating  the  point 
target  problem:  we  are  able  to  define  the  geometry  of  the  problem  in  a  way 
that  allows  relatively  easy  transition  to  the  extended  surface  analysis,  and 
we  can  formulate  analysis  procedures  to  determine  the  form  of  the  "optimum" 
signal  processor.  These  procedures  are  similar  in  the  extended  surface 
analysis. 

The  problem's  geometry  is  defined  in  figure  1.  The  FM/cw  system  is 
instantaneously  located  at  an  altitude  h  above  a  flat  surface  which  is  the 
xy  plane.  The  system  is  moving  with  speed  V  along  a  trajectory  with  an  angle 
of  incidence  6y,  defined  as  the  angle  between  the  flight  path  and  the  vertical 
(z  axis).  At  an  appropriate  later  time  the  system  would  impact  at  point  I. 

The  point  scatterer,  or  target,  is  located  at  point  S  in  the  xy  plane.  Target 
location  is  defined  by  its  position  (ri,  01,  $i)  in  spherical  coordinates 
located  with  origin  at  the  FM/cw  system  at  point  R.  For  the  purposes  of  this 
section,  1  has  only  the  single  value  unity.  However,  we  use  the  notation  so 
that  additional  scatterers  may  be  introduced  in  the  next  section. 

We  define  the  complex  transmitted  signal  of  the  FM/cw  system  as 
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Figure  1.  Geometry  for  FM/cw  system  located  at 
point  R  and  a  scatterer  at  point  S. 
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where 


vT(t)  $  vTe3“Ot+JB<t)  , 

VT  ■  peak  signal  amplitude  , 

i  ^  /T  , 

too  ■  transmitted  carrier  angular  frequency. 


8(t)  =  D  /  m(C)d5  , 


D  ■  modulator's  constant  (rad/s-V) , 
m(t)  ■  modulation  message. 

He  also  make  the  following  additional  definitions: 

&r(°i>  +i)  “  transmit  antenna  one-way  voltage  pattern 
in  direction  (8^,  , 


Gr(9i»  <H)  “  receive  antenna  one-way  voltage  pattern 

in  direction  (8*,  (9 

a l  ”  cross  section  of  scatterer  at  point  S  .  (10 

The  system  of  interest  is  assumed  to  have  the  form  shown  in  figure  2 
where  the  indicated  waveforms  are  the  complex  ones  for  analysis.  The  actual 
system  waveforms  are  the  real  parts  of  the  complex  signals.  The  reference 
signal  that  serves  as  the  local  oscillator  for  the  receive-path  mixer  is 
assumed  to  have  the  arbitrary  form 

vref(t)  -  Z(t)eJ“ot  ,  (11 

where  Z(t)  is  an  arbitrary  complex  envelope 

Z(t)  £  A(t)ejB(t)  .  (12 

We  seek  to  find  A(t)  and  B(t)  so  that  the  response  vjp(t)  is  optimum  in  an 
appropriate  sense. 


The  system  IF  output  can  be  shown  to  be 

ft  It 

Vj\Gt(®£»  ^i)^x(®i»  ^1)^1  j2W(jr£/c+jB(t) 

VH(t>  - - - - A<t)e 


-3b(‘-41) 


,  (13 


where  c  is  the  speed  of  light  in  free  space. 
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Figure  2.  Form  of  FM/cw  system  that  uses 
an  arbitrary  reference  signal. 


(14) 


Next  we  assume  that,  over  the  signal  processing  time  of  interest,  r^ 
can  be  approximated  by  the  first  two  terms  in  its  series  expansion: 


where 


*i(t)  s  roi  ~  *1*  , 


and 


Thus , 


where 


r0^  ■  constant  (value  of  r*  at  t  ■  0) 

(15) 

.  A  dri(t) 

(16) 

r±  dt 

t-o 

2w0ri/c  =  (2o)0r0i/c)  -  (2o>0ri/c)t 

(17) 

A 

”  “0Toi“  udlt  * 

(18) 

T0i  "  2r0i/c 

(19) 

wui  ■  2u0ri/c  . 

(20) 

If  we  now  assume  that  target  motion  affects  mainly  the  return  signal's 
phase  but  has  little  affect  on  the  complex  envelope  (low-frequency  modulations) 
then  (13)  becomes 


vIF(t) 


VV9i»  ♦1)G*(81> 

(4ir)3/2ri2 


j<DOToi->jU)dit+jB(t)-je(t-rOi) 

A(t)e 


(21) 


The  principal  purpose  of  our  analysis  shall  be  to  find  the  power  spectrum 
of  the  IF  output.  In  an  ideal  case  this  spectrum  should  contain  a  single  pair 
of  impulses  at  the  Doppler  angular  frequencies  ±  wdi.  The  power  spectrum  is 
computed  by  Fourier  transforming  the  time-average  of  the  IF  signal's  auto¬ 
correlation  function.  The  autocorrelation  function  of  the  real  part  of  (21) 
is 


RlF(t,  t  +  e) 


V^4(ei,  ♦i)G^(0i,  <fri)*2ai 

2(470  3rj.4 


A(t)A(t  + 


cos  [-wdie  +  B(t  +  e)  -  B(t)  -  B(t  +  e  -  Xoi)  +  g(t  -  Tot)] 

+  cos  [2u»0Toi  -  2udit  -  u)di c  +  B(t)  +  B(t  +  e)  (22) 

-  B(t  -  xoi)  -  8(t  +  e  -  xoi)]l  . 


In  writing  (22)  the  antenna  patterns  have  been  assumed  to  be  real  functions 
No  generality  is  lost  in  this  assumption  because  the  only  change  in  (22) 
for  the  general  case  is  the  subtraction  of  twice  the  phase  angle  of 
<»T(0i,  ♦i)GR(0i,  4i)  from  the  argument  of  the  second  cosine  term,  and,  as  will 
result,  the  second  term  is  nearly  zero  when  the  time  average  is  formed  and  the 
optimum  reference  is  found. 


Ue  form  the  time  average  of  (22)  for  periodic  modulations  and  note  that 
the  second  term  will  nearly  always  be  small  because  of  rapid  cycling  of  the 
cosine  factor  due  to  the  term  2u^^t  in  its  argument.  Ue  therefore  get: 


RifCO  *  A 


j  V  2G  2(0i,  ♦i)GR2(01,  ♦i)X2o1 

*  +  *)>-  - 3 — u - Re 

'  2(4ir)  t± 


I 


T/2 

1  f  A(t)A(t  +  e)eJB(t+e)_je(^_Toi)~jB(t)+JB(t~Toi)dte_:l“die 
T  -T/2 


(23) 


Now  for  a  given  tqi  and  prescribed  0(t),  the  Fourier  transform  of  (23)  will 
possess  impulses  at  ±  only  if  the  integral  is  Independent  of  e .  Furthermore , 
the  "strength"  of  the  impulses  will  be  largest  if  the  integral  is  maximized. 
Since 


b  b 

f  f (x)dx  <  /  | f (x) | dx  ,  (24) 

a  a 

the  integral  (23)  is  maximized  by  first  selecting  B(t)  so  that  the  exponent 
is  zero  for  all  e: 

B(t)  -  B(t  -  T0i)  .  (25) 


The  integral  is  then  further  maximized  by  choosing 

A(t)  -  constant  .  (26) 

The  "optimum"  reference  for  the  system  to  use  for  a  point  target  at  range 
delay  tqi  is  therefore 


Z(t)  -  A(t)ejB(t)  -  eJ8(t"Toi) 


(27) 


where  we  select  A(t)  -  constant  -  1  for  convenience. 


i  ■iff'1*  IT"*1- 

V 


With  the  optimum  Z(t)  we  have 


RlF<e) 


VT2GT2Cflj,  4>i)GR2(9lt  ♦j)X2ff1 
2(4ir)3ri‘* 

VT2gT2(01’  *i>GR2(ei»  4>±)  A2o±Tr 
2<4*)3r14 


cos  (udie) 


[6(w  -  &»di)  +  6(u  +  Wdi)]  •  (29) 


Thus,  the  power  spectrum  possesses  the  desired  impulses  at  the  Doppler  angular 
frequency. 

3.  DISTRIBUTED  FLAT-SURFACE  TARGET 

To  develop  the  distributed  target  problem  we  shall  first  represent  the 
surface  as  a  finite  collection  of  discrete  scatterers  having  random,  statistically 
independent  phase  angles,  uniformly  distributed  on  (0,  2ir) .  The  number  of 
scatterers  is  then  allowed  to  become  infinite  to  form  a  continuous  surface. 
Summations  present  in  the  autocorrelation  function  of  the  IF  output  signal 
become  integrals. 

For  a  finite  number  of  scatterers  the  IF  signal  now  becomes  a  sum  of  terms 
like  (21)  provided  scatterer  cross  sections  are  properly  interpreted.  To  make 
the  subsequent  transition  to  a  continuum  of  scatterers  direct,  we  make  the 
following  substitution: 

/5i  /doi  e^Yi  ,  (30) 

where  doi  is  the  cross  section  of  the  small  (differential)  scatterer  i  that  has 
the  random  phase  angle  Yi.  Cross  section  doi  will  equal  the  product  of  Oo>  the 
cross  section  per  unit  area  that  is  typical  of  the  surface  type,  and  the 
scatterer' s  differential  scattering  area,  dAi,  in  the  xy  plane: 

do*  *  <Tq4Ai  ■  Oori2  tan  8i  dOid^i  .  (31) 

In  general,  Oq  is  a  function  of  the  angle  of  the  incident  wave  that  causes 
scattering.  Thus, 

oq  -  ffo(9i)  ,  (32) 

a  function  of  Q±. 

By  substituting  (30,  (31)  and  (32)  into  (21)  the  IF  output  becomes 

-  Z  VGT(et»  A^^-jYi-ja)dit+JB(t)-3B(t-T01) 

.1: 


The  autocorrelation  function  of  the  real  part  of  (33)  is  found  to  be 


Rip(t,  t  +  e)  -  E[Re  {vip(t)}  Re  {vip(t  +  e)>] 

V  2X2  Gj(6i,  +i)(4(0i,  H) 

-  ,7T-ry  Re  Z  — - r - A(t)A(t  +  t)da± 

K  J  all  i  r±H 

. 


(34) 


#  e-jB(t)+jB(t+e)-ju)<jiC-j$(t+e-Toi)+j0(t-TOi) 


where  E[*]  stands  for  the  statistical  expectation  operation.  Next,  we  drop 
subscripts  i  and  allow  the  sum  to  become  an  Integral  using  (31) .  The  integral 
involves  only  angles  8  and  $,  however,  because  scattering  patch  distance  r 
must  be  restricted  such  that  the  scatterer  lies  in  the  xy  plane.  This  location 
is  preserved  by  substituting 


r  ■  h/cos  0 


(35) 


2r  _  2h 
c  c  cos  8  * 


(36) 


RIP(t,  t  ♦  «)  -  »■  +  e)e‘JB «>**(*•> 

2(4ir)nZ  ( 

w/2  2ir  9 

•  /  Oq (8)  sin  0  cos  0  /  ^(0,  ♦)GR(0,  $)e  ^dx881110008*^  (37) 

6-  0  ♦-« 


.  e  j“dh«:CO80'"jB(t+e-To)+j0(t-To)de 

where  the  radial  Doppler  angular  frequency,  of  any  differential  scattering 
patch  is  related  to  its  components,  u><jh  and  udx>  in  the  h  and  x  directions, 
respectively  (see  figure  1),  by 


with 


ud  “  udh  C08  0  +  udx  8  in  0  coe  $  , 

(38) 

«  (2«0V/c)  cos  0y 

(39) 

*  (2<U(jV/c)  sin  0y  . 

(40) 
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3. 1  Selection  of  Reference  Signal  for  Periodic  Modulations 

Our  ultimate  goal  is  to  find  the  power  spectrum  of  the  IF  output  signal; 
it  is  the  Fourier  transform  of  the  time  average  of  the  autocorrelation  function 
(37).  As  in  the  point  target  case  we  seek  to  make  this  time  average  as  indepen* 
dent  of  e  and  as  large  as  possible.  The  only  factor  of  (37)  involved  in  the 
time  average  is 


A(t)A(t  +  e)e-jB(t)+J0(t-TO)+jB(t+e)-jS(t+e-TO)| 

i  T/2  ' 

/  A(t)A(t  +  s)e^B<t)+^(t-T0)+JB(t+e>-^(t+e-T0)dt  . 

-T/2 

This  integral  will  be  independent  of  e  and  maximum  when 
B(t)  =  8(t  -  t0) 

A(t)  =  constant. 


(41) 


(42) 

(43) 


We  see  that  the  optimum  holds  for  only  one  value  of  tq  [from  (36)  this  means 
only  one  value  of  0],  With  the  optimum  A(t)  and  B(t)  (41)  becomes 


Time  Average 


1  T/2  ■j4"ccos6)+j6(t‘ccos0) 

"  T  -T/2  6 


(44) 


t+e- 


2h 


♦  e 


C  COS  0  0 


,V«(e*n^re) 


dt 


where  0q  is  the  angle  at  which  optimization  occurs. 


To  reduce  (44)  further  we  consider  periodic  modulations  with  fundamental 
angular  frequency  wj.  Now  because  the  modulation  m(t)  is  periodic  by  assumption, 
the  quantity  exp{j0(t  +  e)  -  jB(t)}  is  also  periodic  with  the  same  period  and 
will  have  a  Fourier  series  defined  by 


jB(t+e)-jB(t)  » 


S  Cn(e)e 


jnc^pt 


(45) 


n=-  “ 


where 


(46) 
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On  substitution  of  (46)  Into  (44)  a  simplified  for*  is  found  for  the 
required  time  average,  which,  when  used  in  (37) ,  allows  the  time-averaged 
autocorrelation  function,  denoted  by  Rip(e),  to  be  written  in  the  form 


It  appears  unlikely  that  (47)  can  be  solved  analytically  for  arbitrary 
patterns  and  arbitrary  surface  scattering  function  O(j(0).  Therefore,  in  the 
remainder  of  this  report  we  seek  to  simplify  (47)  in  order  to  obtain  tract¬ 
able  results. 

3.2  Constant  Antenna  Patterns 


The  two  integrals  in  (47)  are  intimately  coupled  through  the  transmit  and 
receive  pattern  functions.  We  may  decouple  the  integrals  to  a  great  extent  by 
assuming  the  patterns  are  constants  over  the  angles  of  interest.  Thus, 

Gt(0,  ♦)  3  Gxo  (48) 
Gr(0»  ♦)  =  Gro  . 


For  many  FM/cw  systems  these  assumptions  are  reasonable. 


By  substituting  (48)  and  (49)  into  (47)  and  making  use  of  a  known 
integral1  we  may  obtain 


„2.2„  2_  2  . 

V-A  Gt  G_  it  *  2  tt/2 

Rxp(e)  -  - — —  2  |Cn(e)|  J  °o(®)  cos  e  8in  e 

(4ir)  ri  n»-  00  0-0 


•  Jo(“dxe  8in  0)  cos  |wdhE  C08  0  “  nuT  T‘tai" 


1 


COS  0n 


:l( 


dO 


(50) 


The  computation  of  (50)  and  the  subsequent  Fourier  transformation  to  obtain 
the  IF  power  spectrum  can  be  done  by  computer.  However,  even  for  relatively 
simple  functions  that  represent  Oq(0),  (50)  is  difficult  to  solve  analytically 


1 

M.  Abramowltz  and  I.  A.  Stegun  (editors),  Handbook  of  Mathematical  Functions 
with  Formulas,  Graphs,  and  Mathematical  Tables  (p  360,  #9.1.21). 
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in  a  manner  that  leads  to  readily  interpreted  expressions . T  In  order  to  obtain 
at  least  one  analytical  solution  to  (50)  we  shall  make  one  additional  simpli¬ 
fying  assumption. 

4.  VERTICAL  INCIDENCE 

If  the  FM/cw  system  is  vertically  incident  on  the  flat  surface,  0y  “  0, 
and  therefore  u»dx  ■  0.  The  autocorrelation  function  (50)  can  now  be  reduced 
further.  If  we  define  functions  Scn(b>)  and  In(ti>)  by 


SCn(a.)  &  Fz{ |Cn(e)j2} 


V“) 


"4 


!  a o  ( 6)  cos  6  sin  6  cos  (u^he  cos  6 


-  “IT (eST  -  cosT^J  M)  ’ 


where  F£{  • }  represents  taking  the  Fourier  transform  with  respect  to  e ,  the  power 
spectrum,  denoted  by  Sjjr(u) ,  which  is  the  Fourier  transform  of  (50) ,  becomes 

vtx^togrow  "  1  " 

SIF(a0  -  fe{RiF(g)}  ■  °  f—  E  /  8-  (5)In(w  -  5)dC  .  (53) 

(Air)  h  n«-«  -  “ 

In  obtaining  (53)  use  is  made  of  the  fact  that  the  transform  of  the  product 
of  two  ,,tlme,,  functions  is  l/2ir  times  the  convolution  of  their  two  spectral  • 

For  our  assumed  problem  where  m(t)  is  periodic  we  have  |cn(e)|2  periodic  in 
e  so  that  it  has  a  Fourier  series  defined  by 

|Cn(e)l2  -  l  Ekn.)k”Ic  ,  (54) 

k—  » 


where 


i  *7* 

Efa-f  t  . 


Peyton  Z.  Peebles,  Jr.,  Probability.  Random  Variables,  and  Random  Signal 
Principles.  McGraw-Hill,  1980  (p  249). 

*  Solutions  in  terms  of  multiple  infinite  series  expansions  can  be  obtained  but 
interpretation  of  results  is  difficult. 


(56) 


The  Fourier  transform,  or  spectrum,  of  |Cn(e)|  becomes 

sCn<“)  ■  ^e(|cn(e)|2J  “  2n  I  Ekn6(u»  -  ku^)  . 

k— 

On  substituting  (56)  into  (53)  we  have 

2  2  2  2 
V  GT0GROir  " 

SIF(o>)  =  — -  °  ■  Z  Z  EknIn(a»  -  kw^)  .  (57) 

(4ir)3h2  n*=-“  k*-» 

Clearly,  SiF(aj)  is  made  up  of  shifted  replicas  of  In(u)  centered  at  multiples 
of  u>F.  If  we  assume 

l“dhl  S  «t/2  ,  (58) 

these  replicas  do  not  overlap  in  frequency.  Furthermore,  we  are  mainly  in¬ 
terested  in  the  frequency  region  (oj(  s  a»j/2  so  that  only  the  term  for  k  *  0 
in  (57)  is  most  important;  if  we  denote  by  the  (Doppler  region)  part 

of  5j;F(u>)  in  the  region  |u|  <  tox/2,  then  for  |u|  <  a»<jh»  it  becomes 

2  2  2  2  / 

W“>  -  +  [v.w  +  '  •  (59) 

I 

Now  from  (46' 

C_n(e)  -  e"jnU>TeCti(-e)  ,  (60) 


|C-n(e)|2  "  |cn(-e)|2 


-n)  1 2  "  I  E0n  1 2 


from  (55).  It  only  remains  to  evaluate  In(k>)  and  I_n(w)  to  reduce  (59) 
to  a  final  form. 

It  is  shown  in  the  appendix  that  (52)  will  reduce  to 


Peyton  Z.  Peebles,  Jr.,  Probability.  Random  Variables,  and  Random  Signal 
Principles .  McGraw-Hill,  1980  (p  256). 


V 


!„(*>)  “  /  o0  (6)  cos  6  sin  el  we"  c  Lcos  8"  cos  80J 


V 

jnw2jir_i - 1 _ I 

+  ne  CLCO80  co8e0J6(u  + 


6(u  -  (D^  cos  0) 


6(u  +  Ujjjj  cos  8)/d6 


C  LW  cose0Jt  o  <  0,  <  „dh 

W°°  [C°S‘‘(^)]  ie3"^~[ra"^]  0 


(63) 


“dh  <  M 


Finally,  we  may  use  (63)  and  (62)  in  (59)  to  obtain  the  IF  output 
signal's  power  spectrum  in  the  Doppler  region: 


(65) 


Several  observations  can  be  made  relative  to  (64).  First,  the  factor 
J 0) |  will  be  maximum  for  ju|  -  w<jh  because  Su^(ta)  must  be  zero  for  |m|  >  u^h* 
Second,  it  is  typical  that  the  scattering  function  a q(')  is  largest  for 
vertical  incidence  when  its  argument  is  zero,  the  case  here  where  |u|  -  udh* 
Therefore,  the  factor  Oolcos_1(|w|  /<i>dh)]  also  helps  maintain  SW{ |(«)  maximum 
at  |u|  -  a)dh.  Furthermore,  Oq  [•  ]  is  a  dominant  term  in  controlling  the  shape 
of  Su>d (**)•  The  faster  og(*)  decreases  with  its  argument,  the  faster  Sif(w) 
decreases  from  its  maximum  for  |m|  smaller  than  udh*  Third,  we  note  that  all 
Egn  *  0  so  that  all  the  cosine  terms  in  (64)  add  in  phase,  producing  a 
maximum  value,  when  |o>|  *  u>dh  If  we  select  cos  8q  ■  1,  or  80  -  0.  Be¬ 
havior  of  the  sum  of  cosines  will  depend  on  the  choice  of  modulation  function 
through  Cn(c)  and  the  sum's  value  is  difficult  to  estimate.  However,  as  I  ml 


>1 


decreases  below  uidh  the  terns  Involving  higher  values  of  n  begin  to  cycle 
rapidly  so  that  SU{j(uj)  is  expected  to  have  rapid  ripples  in  its  amplitude 
as  |ui|  approaches  zero. 

4.1  Exponential  Scattering  Surfaces 

As  an  example  to  illustrate  the  behavior  of  (64)  assume  ao(8)  is  given  by 

Oq(9)  -  ae“btan0t  (66) 

where  a  >  0  and  b  >  0  are  constants.  This  function  is  not  too  unrealistic 
because  it  is  a  reasonable  fit  to  published  data.  For  b  *  1.4  it  fits  data  for 
terrain  forested  by  trees  about  50-ft  high  when  measured  at  X  band.3  Other 
data,1*  for  sea  surfaces  with  winds  from  10  to  20  knots  when  measured  at  fre¬ 
quencies  near  L  band,  are  approximated  by  (66)  when  b  *  8.2. 

With  0q  ■  0,  (66)  becomes 


a0  co8_1[(|u|  /(i>dh)]  -  ae  ' 


while  (64)  reduces  to 


j;  2  2  2 

-  f  ,  V  gtogro*  a  ,  , 

5  (w)  - - - -  |w|e 

d  (4ir)3h2udh^ 


+2  I  Eon  cos 
n*l 


-bVVfv 

|»|e  <E00 

[-“(B  -  ■)]}• 


I“l  <  “dh  » 


l“l  >  “dh  * 


The  precise  values  of  (68)  depend  on  the  coefficients  Eon  which,  in  turn, 
depend  on  the  modulation  employed.  Because  of  this  fact,  only  the  general 
behavior  of  (co)  can  be  estimated.  Based  on  earlier  comments  pertaining 
to  (64)  a  rough  sketch  of  how  (68)  might  behave  is  shown  in  figure  3. 

4.2  Sawtooth  Modulation 

Additional  evaluation  of  (68)  requires  the  modulation  be  specified.  For 
a  sawtooth  FM  that  sweeps  a  band  B(Hz)  in  T  seconds  and  instantaneously  repeats 
sweeps  repetitively,  it  can  be  shown  that 


^M.  I.  Skolnik.  Introduction  to  Radar  Systems .  McGraw-Hill,  1962  (p  525). 

M.  I.  Skolnik,  Introduction  to  Radar  Systems.  McGraw-Hill,  2nd  Edition, 
1980  (p  475). 
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Figure  3.  Rough  plot  of  5Ud(u)  when  the  FM/cw  system 
has  constant  antenna  patterns  and  vertical 
incidence  over  a  flat  surface  with  expon¬ 
ential  scattering  characteristic. 


CU).  (5^1)^- 0,(^1 


jir[e+n(e/BT)-n] 

e 


(69) 


+  (e/BT)Sa{ir(e  -  n  -  BT)e/BT)eJir^e+n(e/BT)"BT^ ,  0  <  e  <  T  , 


where 


e  »  Be  ,  (70) 

and  we  define  the  sampling  function  Sa(*)  by 

Sa(£)  -  Cl  sin  «)  .  (71) 

Thus, 

+l4f\2Sa2{it(t  -  n  -  BT)e/BT) 

\BT/  (72) 

+  2(gyBFl)(M)Sa|1f(g  "  n)(BTBT  ^|Sa{ir(g  “  n  “  BT)e/BT} 

•  cos  [w(BT  -  n)]  •  0  <  e  <  T  . 


The  first  term  in  (72)  has  Its  largest  amplitudes  when  e  is  small;  it  peaks 
to  largest  values  when  E  -  n.  If  BT  is  reasonably  large  (say  20  or  more)  then 
the  first  term  will  be  fairly  small  for  e  near  BT.  The  second  term  has  its 
largest  amplitude  coefficient  (e/BT)2  when  e  is  near  BT.  However,  since  the 
argument  of  the  sampling  function  is  always  negative  (zero  only  for  E  ■  BT  and 
n  ■  0)  for  n  >  1 ,  the  second  term's  total  amplitude  will  tend  to  be  small  for 
all  E  when  n  i  1.  This  fact  also  means  the  third  term  in  (72)  is  small. 

Thus,  for  large  BT» 


M) 

= ° 


(73) 


<  E  <  BT  . 


On  substitution  of  (73)  into  (65)  we  have 


0.0  «-r  OTW  V  vlv".  v’a'a  a  ’.V-  a 


fBT  -  n)  ’(M-o)2/®!  ,  »T  _  n 

"to  2  '  «t  Sa  Mix  *  "TbtF  ■  lini 


-nir(BT-n)  /BT 


s  0  , 


BT  -  1  <  n 


The  power  spectrum  now  becomes 


2  2  2  2  - 
V„A  G_0„  w2a 


“dh\2 


(4ir)  h  todh  BT 


BT-1 

1  +  2  £  (1 


n  \  f  2h/“dh  ^ 

It)  cos  prTyJf  ‘  1 


It  can  be  shown  that 


N 

2  £  (1  -  cos  (nx)  ■ 

n-1  '  ' 


,  .  i  fain  (Nx/2)12 
N  sin  (x/2)J 


and  (.75)  reduces  to 


sit 


2  2  2  2  2 
V  Gt,gr.'  «  ImI 
(4«)V.-,  “ah 


.llt2  BT“fh  ( “ah 
*ln  .  c  \1U 


(BT)2  sin2 


‘^/“dh  ,V 

.« lw ' ) 


:(2“ah) 


where  we  define 


rect  {ifi)  "  l»  l“l  <  w 

-  0,  |u|  >  W  , 

as  usual.  Nov  In  general,  if  we  define 


■ 

*•  I 


then  for  most  systems  x^  <  <  T  and  the  argument  of  the  denominator  sine  term 
in  (77)  is  small  for  all  |u>|  near  Finally,  we  have  the  approximation 


2  2  2  2  2 
V  SAo; a 
aioV^ 


{*BT»(p  -  *)}  ’ 


(2“dh) 


The  normalized  value  of  (80)  is  plotted  in  figure  4  for  Bx^  *  2.  The 
constant  terrain  curve  presumes  b  ■  0  while  the  other  curve  presumes  b  ■  1.4. 
For  comparison,  figure  5  shows  similar  plots  for  Bx^  -  4.  Clearly,  the  power 
spectrum  tends  to  spread  as  the  system  becomes  closer  to  the  scattering 
surface,  for  a  given  bandwidth  B. 


Brh  -  4 

(CONSTANT 

TERRAIN 


eV. V*" . tij^^r-^'VwV'iV  •  . 


APPENDIX 


The  steps  leading  from  (52)  to  (63)  are  outlined  in  this  appendix. 
From  well-known  transforms2  we  first  write  (52)  in  the  form 


cos  0°J5(u)  -  a)dh  cos  6)d6 


«/2  jnuj  — - i - ±— -I 

+  /  oq(0)  cos  0  sin  6 ire  C  ^cos  0  cos  ®0-l6((u  +  u>dh  cos  0)d0  . 

0 

Since  0  can  have  values  only  from  0  to  rr/2  the  impulse  function  in  the  first 
integral  above  will  "occur"  only  for  0  <  u  <  Wdh*  More  precisely,  it  occurs 
when 


!»<“) 


ir/ 

I 


/•% 


a0(6)  cos  0  sin  0  ire 


cos  0 


0  -  01  -  cos”1  (u/ajjh) . 

For  the  values  of  <■>  for  which  the  first  term's  impulse  can  occur,  the  second 
term's  impulse  cannot  "occur"  leading  to  the  second  integral  above  being  zero 
for  0  <  o>  <  udh*  Similarly,  the  second  term  is  nonzero  for  -t#dh  <  w  <  0  while 
the  first  Integral  term  is  zero.  When  |u|  >  u>dh  neither  impulses  "occur"  for 
any  0  <  0  <  ir/2  and  both  Integrals  are  zero.  Hence  In(w)  *  0  for  |u>|  >  «dh* 

Consider  the  evaluation  of  the  first  Integral  above,  which  is  nonzero  only 
for  0  <  o)  <  Udh*  Define 

f(0)  -  u  -  u,jh  cos  0 

f '  (0)  -  “  “dh  sin  0 

|f'(0)|  -  |u>dhl  lain  0j  -  t>>dh  sin  0 
for  0  <  0  <  ir/2.  We  apply  a  theorem  of  Friedman:5 


2 

Peebles,  Peyton  Z,  Jr.,  Probability.  Random  Variables,  and  Random  Signal 
Principles.  McGraw-Hill ,  1980  (p.  256). 

Friedman,  B. ,  Principles  and  Techniques  of  Applied  Mathematics,  John  Wiley, 
1956  (Fifth  printing ,"1962)  pp.  lSS-TIT 
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VIV, 


Theorem  If  $(x)  is  continuous  at  x  ■  x<j 
and  f(x)  is  a  monotonic  function  of  x 
that  vanishes  at  x  ■  xq ,  then 

*  *(x0) 

t  *(x)«[f(i)]dx  . 

We  identify  8  with  x,  81  with  x0 ,  and  apply  the  theorem  to  the  first  integral 
term  in  In(<a)  by  letting 

....  4  ...  J  - 1^9?] l 

♦(8)  *  <*o(®)  cos  6  sin  ®ve  u  J  C  • 


Thus, 


I„Ca.)  -  ^2 

“dh 


J  \ 

cos  0  sin  0<ire  u  J  l  * 

,  2hfudh  1  1 

T  -1/  »  “JnUiTc  L  «  "  cos  0O  J 

On  I  COS  I - lie  u  » 

0 1  Wh/J 


for  0  <  id  <  wdh* 

By  repeating  the  above  procedure  by  letting  u>  *  -|oa|  for  -Wdh  <  w  <  0 
and  identifying 


.  .  .  fl{  ^XUl*^l|coe0-  cos  0(J  } 

♦(0)  ■  oo(0)  cos  0  sin  0jxe  u  J  I 

f(0)  ■  -|ui|  +  ujk  cos  0,  0  <  ]u|  <  a>dh 

01  -  cos_1(|«|/wdh) 

f’(0l)  ■  ~“dh  8ln  ®1  » 

we  develop  the  second  Integral  for  I^Cu)  as 

InM  -  ^  , 


for  -w^h  <  w  <  0. 


By  combining  the  above  two  results  for  In(u)  we  obtain  (63) . 
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